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case of Caratheodory Systems and Networks 



- 1—1 

X 



Mario di Bernardo ' 6 Davide Liuzza a Giovanni Russo a 

. a Department of Systems and Computer Engineering, 

CN \ University of Naples Federico II, 80125, Naples, Italy. 

■ {mario . dibernardo , davide . liuzza, giovanni . russo2}@unina. it 

b Department of Engineering Mathematics, University of Bristol, BS8 1TR, 
^-j. ■ Bristol, U.K. {m.dibernardo@bristol.ac.uk} 

u 
o 

■i— > . 

In this paper we extend to a generic class of piecewise smooth dynamical systems a 
fundamental tool for the analysis of convergence of smooth dynamical systems: con- 
traction theory. We focus on switched systems satisfying Caratheodory conditions 
for the existence and unicity of a solution. After generalizing the classical definition 
. of contraction to this class of dynamical systems, we give sufficient conditions for 

ly-) \ global exponential convergence of their trajectories. The theoretical results are then 

applied to solve a set of representative problems including proving global asymp- 
totic stability of switched linear systems, giving conditions for incremental stability 
of piecewise smooth systems, and analyzing the convergence of networked switched 
linear systems. 



Abstract 



H ! 1 Introduction 

C3 



Piecewise-smooth dynamical systems are commonly used in Nonlinear Dy- 
namics and Control to model devices of interest and/or synthesize discontin- 
uous control actions e.g., [Cor08], [dBBCK08]. Despite the large number of 
available results on their well-posedness and stability, there are few papers in 
the literature where the problem of assessing their incremental stability and 
convergence properties is discussed. 

The available results usually refer to specific classes of systems and rely on 
ad hoc continuity assumptions. For example, the problem of proving conver- 
gence for piecewise affine continuous systems and networks is addressed in 
[PvdWN05], [PvdWN06], [PPvdWN07], [PvdW08] using a Lyapunov based 
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approach. The methodology extends the approach of Demidovich for smooth 
dynamical systems expounded in [PPvdWN04] . 

Contraction theory is an alternative approach used to study convergence be- 
tween trajectories in smooth dynamical systems (see [LS98] and references 
therein). The contraction approach is based on finding some metric under 
which the matrix measure of the system Jacobian can be proved to be definite 
negative over some convex set of phase space of interest. It can be shown that 
Demidovich approach is related to proving contraction using Euclidean norms 
and matrix measures. Nevertheless, using contraction it suffices to find some 
measure to study the Jacobian properties including non- Euclidean ones (e.g., 
//i, /loo etc). 

Contraction theory has been used in a wide range of applications. For example, 
it has been shown that contraction is an extremely useful property to analyze 
coordination problems in networked control systems such as the emergence of 
synchronization or consensus ([LS98,PS07,RdBS10a,WS05,RdB09a,RdB09b,RS10]). 
Indeed, all trajectories of a contracting system can be shown to exponentially 
converge towards each other asymptotically. Therefore as shown in [WS05], 
this property can be effectively exploited to give conditions for the synchro- 
nization of a network of dynamical systems of interest. Recently, it has also 
been shown that the use of non Euclidean matrix measures can be used to con- 
struct an algorithmic approach to prove contraction [RdBSlOa] and to prove 
efficiently convergence in biological networks [RdBSlOb]. 

Historicallyideas closely related to contraction can be traced back to [Har61] 
and even to [Lew49] (see also [PPvdWN04], [Ang02], and e.g. [LS05], [Jou05] 
for a more exhaustive list of related references). For autonomous systems 
and with constant metrics, the basic nonlinear contraction result reduces to 
Krasovskii's theorem ([SL90]) in the continous-time case, and to the contrac- 
tion mapping theorem in the discrete-time case ([LS98], [BT89]). 

Despite the usefulness of contraction theory in applications, there is no con- 
sistent extension of this approach to the large class of piecewise smooth and 
switched dynamical systems. In [LSOO], it is conjectured that, for certain 
classes of piecewise-smooth systems, contraction of each individual mode is 
sufficient to guarantee convergence of all the system trajectories towards each 
other, i.e. contraction of the overall system of interest. Also, in [ERS06], it is 
noted that contraction theory can be extended to a class of hybrid systems 
under certain assumptions on the properties of the reset maps and switching 
signals. 

The aim of this paper is to start addressing systematically the extension of 
contraction theory to generic classes of switched systems. Here, we focus on two 
types of switched systems of relevance in applications: (i) piecewise-smooth 
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continuous (PWSC) systems (a class of state-dependent switched systems), 
and (ii) time-dependent switched systems (TSS). The goal is to obtain a set 
of sufficient conditions guaranteeing global exponential convergence of their 
trajectories. The extension of contraction theory to Filippov vector fields will 
be presented in a separate paper currently in preparation. 

From a methodological viewpoint, to investigate contraction properties of 
PWSC and TSS we focus on systems whose vector fields satisfy Caratheodory 
conditions for the existence and uniqueness of an absolutely continuous solu- 
tion (see [Fil88] and [Cor08] for further details). We prove that, as conjectured 
in [LSOO], for this class of systems, contraction of each individual mode suffices 
to guarantee convergence of all the system trajectories towards each other, i.e. 
contraction of the overall system of interest. Contrary to the results presented 
in [PvdWN05], [PvdWN06], [PPvdWN07], [PvdW08], we do not require find- 
ing incremental Lyapunov functions for the system of interest. We use instead, 
as commonly done for analysing contraction in smooth systems, a generic con- 
dition on the existence of some metric in which the Jacobian of each mode of 
the system of interest is definite negative. We then apply the theoretical results 
to study convergence of some representative examples, including a network of 
time-switched systems. 

A preliminary version of some of the results presented in this paper were pre- 
sented at the 18 th IFAC World Congress in Milan in September 2011 [RdBll]. 



2 Mathematical Preliminaries 



Before presenting the main results of the paper, we introduce here some es- 
sential definitions and notation that will be used in the rest of the paper. 

Let x be an n-dimensional vector. We denote with |x| the norm of the vector. 
Let A be a (real) matrix. Then, \\A\\ denotes the norm of A. We recall (see for 
instance [MLH08]) that, given a vector norm on Euclidean space (|-|), with its 
induced matrix norm \\A\\, the associated matrix measure /i is defined as the 
directional derivative of the matrix norm, that is, 

^ :=S55<P + M||-1). 

For example, if |-| is the standard Euclidean 2-norm, then n{A) is the maximum 
eigenvalue of the symmetric part of A. As we shall see, however, different norms 
will be useful for our applications. Matrix measures are also known as "loga- 
rithmic norms" , a concept independently introduced by Germund Dahlquist 
and Sergei Lozinskii in 1959, [Dah59,Loz58]. The limit is known to exist, and 
the convergence is monotonic, see [Str75,Dah59]. 
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In what follows we report the analytic expression of some matrix measures 
used in the paper: 

• fi^A) = max, (ajj + J2i^j Wijl); 

• pL 2 (A) = max; \i (j [A + A T )); 

• fioo{A) = max; (an + £ j¥i 

More generally, we will also make use of matrix measures induced by weighted 
vector norms, say |x| ei = |Ox|j, with a constant invertible matrix and 
% = {l,2,oo}. Such measures, denoted with fie,i, can be computed by using 
the following property: HQ ti (A) = \ii (0A6 -1 ), \/i = {1, 2, 00}. Obviously, any 
other measure can be used. 

We will also use the following definitions. 

Definition 1 (K-reachable sets [RdBSlOb]) Let K > be any positive 
real number. A subset C C W 1 is i^-reachable if, for any two points x and yo 
in C there is some continuously differentiable curve 7 : [0, 1] — > C such that: 

(1) 7 (0) = a;o ; 

(2) 7 (1) = y and 

(3) \i{r)\<K\y -x \,Vr. 

For convex sets C, we may pick 7(7-) = Xo + r(y — x ), so 7'(r) = y — x and 
we can take K — 1. Thus, convex sets are 1-reachable, and it is easy to show 
that the converse holds as well. 

Definition 2 (Flow of a Dynamical System) Given a dynamical system 

x = f(x,t), xeVCW n ,teR, 

we define its flow ip(s, t , x) : IR + x R+ x R™ — > R" as the operator (see [Son98] 
for a definition) such that 



*o, X) = f(v( s , to, X),t), <p{0, k, x) = X- 

In applications of the theory, it is often the case that V will be a closed set, for 
example delimited by some hyperplane in the phase space, which could e.g. 
model constraints on the state variables of the system. We remark here that 
for a non-open set, differentiability in x means that the vector field f(»,t,) 
can be extended as a differentiable function to some open set which includes 
V, and any continuity hypotheses with respect to (t, x) hold on this open set. 
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2.1 Caratheodory Solutions 



We give the following preliminary definitions: 

Definition 3 A function g(t) : [to, +00] — > M. n is said to be measurable if, for 
any real number a, the set {t G [to, +00] : g{t) > a} is measurable in the sense 
of Lebesgue. 

Definition 4 A function l(t) : [to, +00] — > K. is summable if the Lebesgue 
integral of the absolute value of l{t) exists and is finite. 

Definition 5 A function z(t) : [a, b] — > lR n is absolutely continuous if for all 
e > there exists 5 > such that for each finite collection [a\, bi] . . . [a n , b n ] of 
disjoint sets in [a,b], it holds that J2k \t>k — Q>k\ < $ =^ J2k \ z (bk) — z(a,k)\ < s. 

Now we are able to give the definition of a Caratheodory solution of a differ- 
ential equation (see also [Cor08] and references therein): 

Definition 6 (Caratheodory solutions) Let us consider a domain V C 
W 1 and a dynamical system of the form: 

x(t) = f{x(t),t), x{t ) = x , (1) 

where f : T> x M + — > W l . A Caratheodory solution for this system is an 
absolutely continuous function x(t) that satisfies (1) for almost all t e [to,^i] 
(in the sense of Lesbesgue), where \po,t]\ is an interval where the solution x(t) 
is defined. That is, a Caratheodory solution of (1) is an absolutely continuous 
function x(t) such that: 

x(t) = x(t ) + /* f{x(r),r)dr, t G [t , t ± ]. 

J t 

In the common cases where the solution x(t) can be extended forward in time, 
the same definition holds for each ti > t . 

An useful result provides sufficient conditions for the existence of a Caratheodory 
solution of the system (1). 

Theorem 1 (Existence and Uniqueness of Caratheodory solutions) A 

Caratheodory solution of system (1) exists if: 

(1) for almost all t G [0, 00] ; the function x >->■ f(x,t) is continuous for all 
x G V; 

(2) for each x G V, the function t i-> f(x,t) is measurable in t; 
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(3) for all (x,t) 6Dx [0,+oo] 7 there exist 5 > and a summable function 
m(t) defined on the interval [t,t+ 5] such that \f(x,t)\ < m(t). 

Moreover, the solution is unique, if the following additional condition is sat- 
isfied: 

(4) (x — y) T (f(x, t) — f(y, t)) < l(t)(x — y) T (x — y), where l(t) is a summable 
function. 

Notice that, as discussed in [Fil88], p. 10 and proved in [AB58], equations 
that satisfy the above (Caratheodory) conditions and those required for the 
uniqueness of a solution show continuous dependence on initial conditions. 

2.2 Caratheodory Systems 

Following [dBBCK08] p. 73, we define a piecewise smooth dynamical system 
as follows. 

Definition 7 Given a finite collection of disjoint, open and non empty sets 
Si, . . . S p such that V = Uf=i Si C R n is a connected set, a dynamical system 
x = f(x,t) is called a piecewise smooth dynamical system (PWS) if it is 
defined by a finite set of ODEs 

f(x,t) = Fi(x,t) xeSi. 

The intersection E^- := Si dSj is either a lower dimensional manifold or it is 
the empty set. Each vector field Fi(x,t) is smooth in both the state x and the 
time t for any x e Si. Furthermore it is continuously extended on the boundary 
dSi. 

We now introduce two important classes of switched systems that will be 
analyzed in the paper and for which a Caratheodory solution exists. 
A piecewise smooth dynamical system is said to be continuous (PWSC) if the 
following two conditions hold: 

(1) the function (x,t) h- > f(x,t) is continuous for all x G M. n and for all t > t ; 

(2) the function x >->■ Fi(x,t) is continuously differentiable for all x G Si and 
for all t > t . Furthermore the Jacobians jj£(x,t) can be continuously 
extended on the boundary dSi. 

Notice that in order for condition (1) to be satisfied the functions (x,t) h- > 
Fi(x,t) must be continuous for all t > t and x G S iy and, for all x G E^- ^ 
and all t > t , it must hold Fi(x,t) = Fj(x,t). 
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We now give the definition of a time-dependent switching system according 
to [Lib03]. 

Definition 8 A time-dependent switching system is a dynamical system of 
the form 

x = f(x,t,a), xGl", (2) 

where a(t) : [0, +00) — > £ = {1,2,..., p} is a piecewise continuous time- 
dependent switching signal taking one over p finite possible values. 

Note that according to this definition, we are excluding the case where infinite 
switchings occur over finite time so that Zeno behavior cannot occur (see 
[Lib03] for further details). 

We will refer to Caratheodory systems to indicate any PWS or TSS satisfying 
the conditions for the existence of a Caratheodory solution given in Theorem 1. 



3 Contraction Theory: a brief overview 



In this section, the notion of contraction is briefly summarized for a generic 
nonlinear system of the form: 

x = h(x, t) (3) 

where h : T> x R 1— > IR ra is an m-dimensional vector field assumed to be suffi- 
ciently smooth. (For further details see [LS98,RdBS10a,RdBS10b].) 

Definition 9 The smooth dynamical system (3) is said to be (infinitesimal) 
contracting over a K -reachable set C C W 1 , if there exists some matrix mea- 
sure, ji, and a positive scalar c, such that 



dh . , 



< -c, VxGC,Vt>0. 



The positive scalar c is said to be the contraction rate of the system. 

The basic result of nonlinear contraction analysis states that, if a system is 
contracting, then all of its trajectories exponentially converge towards each 
other, see [LS98,RdBS10b]. 

Theorem 2 (Contraction) Assume that (3) is contracting and let x(t) and 
x(t) be any two of its solutions with initial conditions x(to) = Xo £ C and 
x(t ) = Xq G C. Then, for any t > t , it holds that 

\x(t) - x(t)\ <K\x - x \ e~ ct ^ to) 
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4 Contraction of Caratheodory Systems 



Contraction theory has been properly studied mostly in the case of smooth 
nonlinear vector fields. The case of switched and hybrid systems is only marginally 
addressed in the existing literature [LSOO], [ERS06]. 

In this section, we seek sufficient conditions for the convergence of trajectories 
of PWSC systems (a generic class of systems with state-dependent switch- 
ings) and time-dependent switched systems. For the sake of clarity, we keep 
the derivation for the two cases separate. 



4-1 Contraction of PWSC systems 



We start with PWSC systems as defined in §2.2. We can state the following 
result: 

Theorem 3 Let C C V be a K-reachable set. Consider a generic PWS system 
of the form 



x = f(x,t) = < 



F x (x,t) xeSi 

i (4) 
F p (x, t) x e S p 



defined as in Definition 7 for all x e C and with smooth manifolds for all 
i, j = 1, . . . ,p. Suppose that: 

(1) it fulfills conditions for the existence and uniqueness of a Caratheodory 
solution given in §2.1; 

(2) there exist a unique matrix measure such that 

"(^ (x '* ) )-" Ci ' 

for all x e Si and all t > t , with Ci belonging to a set of positive scalars 
(in what follows, we will define c := mhij Ci). 

Then, for every two solutions x(t) — (p(t — to, to, £) and y(t) — (p(t — to, to, Q 
with £, C G C, it holds that: 

|x(t)-2/(t)|<Ke- c ^|e-C|, 

for all t > t such that x(t), y(t) EC. If C is forward-invariant then all trajec- 
tories rooted in C converge exponentially towards each other. 
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Proof. Given two points x(to) = £ and y(to) = C an d a smooth curve 7 : 
[0, 1] 1 — y C such that 7(0) = £ and 7(1) = (, we can consider tf)(t, r) := 
Lp(t — t ,t , 7(r)) as the solution of (4) rooted in i[)(t , r) = 7(r), with r G [0, 1]. 
Notice that ip(t, r) is continuous with respect to r for all t. Notice also that 7 
can be chosen so that ip(t,r) is differentiable with respect to r for almost all 
the pairs (t, r). Let 



a.e. in t, a.e. in r. 



(5) 



Thus we have: 




^;(fm,r),t)), 



a.e. in t, a.e. in r, 



In what follows we will use the shorthand notation a. e. to denote the validity of 
a given expression almost everywhere in both t and r, unless stated otherwise. 

Since 

we can write: 

d 

—w(t,r)=A(?p(t,r),t)w(t,r), a.e., (6) 

where we have denoted by A(x, t) the Jacobian of the PWSC system (4), which 
can be defined as: 



A(x,t) = yt(x,t) 



f±(x,t), VxeSi, 



for almost all the pairs (x, t) apart from those points where x G S^, for some 

The next step is to show that the solution t i->- w(t,r) of (6) is a continuous 
function for any fixed r G [0, 1]. 

Indeed, without loss of generality, consider the image of the curve 7 under the 
action of the flow <p for a time T such that the system trajectory rooted in 
7 has either crossed the boundary once or it has not (in the case there are 
multiple switchings between t and T, the same reasoning can be iterated). 
Furthermore, let us call r(r) G]to,^[ the time instant at which the trajectory 
eventually crosses the boundary. Suppose, without loss of generality, that at 
t = r(r), the flow switches from region S\ to region £2. Then, we have: 
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ij(t,r) 



<pi(t - t ,t ,ip(t ,r)) 



to < t < r(r), 



<p 2 (t - r(r), r(r), Vl (r(r) - t , t , ^(t , r))) r(r) < t < T. 



Now, to show continuity of w(t,r) with respect to time, from (5) we need to 
evaluate the derivative of ip{t,r) over the interval ]to>^[- We have: 



dr [<fi(t-to,to,^(t ,r))} 
Or )9_[ (fi2 ( t - r(r), r(r), ^(r(r) - t , *o, ^(*o, 0))] r(r) <t <T, 



Continuity of w(t, r) is then guaranteed if 



to<t< r(r), 
(7) 



o-[Vi(* - *o, *o, ^(*o, r))] = lim — [v?2(t-to-r(r),r(r),v?i(r(r)-t ,to,^(to,r)))], 

(8) 



We have 



— ipiis to x) = dX 
dr ' ' S"x or 



(9) 



with s := t — t and \ = V^Oi r ) := V* - Hence, taking the limit t — > r(r) , 
the left-hand side of (8) can be written as: 



W ^)-to,t Q ^) — 



(10) 



Also 



where 



d j. dcp 2 ds dtp 2 di dip 2 dx 

-Ms(t,r) M r), X (r)) = —- + + — - 



s(t,r) 



t-r{r) 
r(r) 

(f! (-s(t ,r),t ,ip ) 



(11) 
(12) 

(13) 



Now, we observe that 



d(p 2 
ds 



F 2 (ip 2 (s(t,r),r(r),x(r)),t) 



dx dipt diptdijj 

— = t (r) H 

dr ds dip dr 



(14) 
(15) 
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where 

^- = F 1 (M-s(t ,r),t ,ij°),t) 

Taking the limit for t — > r(r) ± , we have: 

lim s(t, r) = 

t->-r(r)+ 

hence, since 

<Pi(0,T(r),x(r),t) = x(r) (16) 
we then obtain that (14) yields in the limit 



as 



F 2 ( X (r),r(r)) = ^(^(-s^o.r),^, V°),r(r)) 



Moreover, in the same limit, t — > r(r), from (16) we have: 

dt di 

and 

dx dx 

Therefore, the right-hand side of (8) can be written as 
- F 2 (M-s(t ,r),t ,^),r(r))r'(r) 

+ FiM-^o, r), to, A r(r))/(r) + |^ (17) 

From the assumption that the system vector field is continuous when t = 
r(r), continuity of w(t, r) with respect to t is then immediately established by 
comparing (10) and (17). 

Now, we turn again our attention to equation (6). Fixing r to any value be- 
tween and 1, the Jacobian can be calculated and (6) can be solved to obtain 
(in the sense of Lebesgue): 

w (t + h,r)= w(t, r) + J t t+h r),ti)w(&, r)dfi = 

= w(t, r) + A(ip(t, r),t)w(t, r)h+ 

J t t+h {A{^{$,r),tf)w{$,r) - A(iP(t,r),t)w(t,r))d$, a.e. t, 
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with h being a positive scalar. 

Thus, from the above expression we have 



\w(t + h,r)\ < \\I + hA(ip(t,r),t)\\ \w(t,r)\ 

rt+h 

+ J \A(i/}(-&,r),-&)w{-&,r)-A(il;(t,r),t)w(t,r)\d'd (18) 

Then, subtracting \w(t, r)\ from both sides of the equation and dividing by h 
we obtain 

i (\ w ( t + h,r)\ - \w(t,r)\) <l(\\I + hA(i;(t,r),t)\\ - 1) \w(t,r)\ + 
y t t+h \A(i(j($,r),$)w($,r) - A(iP(t,r),t)w(t,r)\d$, a.e. t, 

Thus, taking the limit as h \ yields: 

^\w(t,r)\<-c\w(t,r)\, a.e. 

Notice that the above expression holds for all those pairs t and r where the 
Jacobian A(-) is defined. Let now M(t) := —c(t—t ), from the above expression 
it follows that: 

j t (\w(t,r)\e- M U)<0, a.e. 

Now, since e~ M ^ is an increasing function and since the function t h-> w(t,r) 
is continuous, the above inequality implies that: 

\w(t,r)\ < \w(t ,r)\e- c(t ' to) < K \£ - (\ e'^^. 

As the function ip(t, r) is continuous and, for all t, the function w(t, r) is 
defined for almost all r, we have: 

^(t,l) -ip(t,0) = f 1 w(t,s)ds. 
Jo 

Thus: 

\ x (t)-y(t)\<K\^-C\e-^- t0 \ 
and the theorem remains proved. 

Obviously, if C is forward- invariant, then trajectories rooted in C will expo- 
nentially converge towards each other. 
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4-2 Contraction of TSS 



The conditions used to prove contraction of PWSC systems can be immedi- 
ately extended to generic systems affected by time-dependent switchings as 
detailed below. 

Theorem 4 Consider an invariant K-reachable set C C W 1 and a time- 
dependent switching system as in (2). Suppose that: 

(1) it fulfills conditions for the existence and uniqueness of a Caratheodory 
solution given in $2.1; 

(2) the function (x,t) h-> f(x,t,a) is continuous for all x G C, for all t > t 
and for all a G S; 

(3) the function x >->■ f(x,t,a) is continuously differentiable for all x G W 1 , 
for all t > to and for all a G S; 

(4) there exist a unique matrix measure such that 



for all x G C, for all t > t and for all a G S, with c a belonging to a set 
of real scalars (in what follows, we will define c 2 := minors c 2 a ). 

Then, for every two solutions x(t) = (p(t,t ,£) and y(t) = (p(t,t ,() with 
C G C, it holds that: 



Proof. The proof follows similar steps to that of Theorem 3. In particular, 
given points x(t ) = £ and y(t ) = ( and a smooth curve 7 : [0, 1] (->■ C such 
that 7(0) = £ and 7(1) = (, we can consider ip(t,r) := tp{t — t ,^0)7( r )) as 
the solution of (2) rooted in i(j(t , r) = 7(7"), with r G [0, 1]. Let 




x(t)-y(t)\<Ke-^^-t\, 



w(t, r) : 



a.e. t. 



dr 



As in the proof of Theorem 3, we can write: 



dw 



dt 




a.e. t 



and then 
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— w(t,r) = A(ip(t,r),t,a)w(t,r), a.e. t, 



with A = -^f(x,t,a) being the Jacobian of f(x,t,a) for almost all t > to- 
Notice that, differently from the case of PWSC systems, here the Jacobian is 
discontinuous only with respect to time t due the fact that the function a(t) 
is piecewise constant. However, we can show that the function t h-> w(t,r) 
is continuous by considering again (10) and (17). In this case, the switching 
instant r(r) is independent from r and therefore all terms containing r'(r) in 
(17) cancel out. The equality of (10) and (17) then immediately follows and 
the rest of proof becomes identical to that of Theorem 3. 

Remarks 

• Our results on the contraction of PWSC systems can be interpreted follow- 
ing the approach presented in [CFPT06] where the asymptotic stability of 
piecewise linear systems obtained by the continuous matching of two stable 
linear systems is discussed. Specifically, under the conditions of Theorem 3, 
we can state that the continuous matching of any number of nonlinear con- 
tracting vector fields is also contracting. As in the case of asymptotic stabil- 
ity discussed in [CFPT06], guaranteeing incremental stability of switched 
systems is not trivial, even when they are obtained by continuously match- 
ing contracting vector field. Thus, the sufficient conditions derived in this 
paper can be useful for the analysis of incremental stability in switched 
systems and the design of stabilizing switched control inputs. 

• The results reported so far do not include the case of differential equations 
with discontinuous right-hand side or Filippov systems where so-called slid- 
ing motion is possible [dBBCK08]. This is ongoing work that will be reported 
elsewhere [dBLR]. 



5 Partial Contraction of PWSC and TSS 

In the previous section, contraction theory was extended to piecewise continu- 
ous and time-dependent systems. In both cases, we were able to show that the 
contraction still implies asymptotic convergence between trajectories. Often 
in applications, it is desirable to prove (or certify) that, at steady state, all 
trajectories of a given system exhibit some property regardless of their initial 
conditions. 

In the case of smooth dynamical systems, the concept of partial contraction 
was introduced in [WS05] to solve this problem. The idea is to introduce an 
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appropriately constructed auxiliary or virtual system, embedding the solutions 
of the system of interest as its particular solutions. If the virtual system is 
proved to be contracting, then all of its solutions will converge towards a 
unique trajectory. In turn, this imples that all trajectories of the system of 
interest, embedded in the virtual system by construction, will also converge 
towards this solution. 

The most notable application of partial contration is its use to prove conver- 
gence of trajectories of all nodes of a network of dynamical systems towards 
each other as for example is required in synchronization or coordination prob- 
lems. In that case, the virtual system is constructed so that trajectories of the 
network nodes are its particular solutions. Proving contraction of the virtual 
system then implies convergence of all node trajectories towards the same syn- 
chronous evolution (see [WS05], [SWR04], [RS10], [RdB09a], [RdBSlOa] for 
further details and applications). 

Using the extension of contraction to switched Caratheodory systems pre- 
sented in this paper, we can also extend partial contraction to this class of 
systems. For example, we can prove the following result for PWS systems. 

Theorem 5 Consider system a PWS of the form (4) and assume that there 
exists some system 

y = v(y,x,t) (19) 

such that: 

• v(x,x,t) = f(x,t); 

• v(y, x, t) is contracting in the Caratheodory sense with respect to y and for 
any x. 

Let y be the unique solution towards which all trajectories of (19) converge. 
That is, there exists some c > such that, for any solution of (19), say y(t): 

\y(t)-y(t)\<\y(to)-y(t )\e- c{t - to) 

Then, all the solutions of (4) converge towards y, i.e. 

\x(t) - y(t)\ ->■ 0, t ->■ +oo 

System (4) is said to be partially contracting while system (19) is termed as 
virtual system. 

Proof. Indeed, we only need to observe that by construction any solution 
of (4), say x(t), is also a solution of the virtual system. Now, since (19) is 
contracting, then all of its solutions will converge towards y. This in turn 
implies that 

\x(t) - y(t)\ ->■ a.e. 
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as t — > +00. □ 



We remark here that Theorem 5 allows to prove that all the solutions of the 
real system of interest converge, at steady state, towards a unique solution 
even if it is not contracting. The key point of such a result is indeed that 
of constructing a contracting system which embeds the solutions of the real 
system. The main implication of this is that trajectories of the real system will 
converge towards each other but convergence will not be exponential and, in 
general, it may be non-uniform. That is, Theorem 5 only ensures that, after a 
sufficiently long time, distances between all the solutions of (4) shrink. In some 
special case, see e.g. [RS10], the dimensionality of the virtual system is lower 
than that of the real system of interest: this is typically the case of systems 
with symmetries, such as Quorum Sensing networks. A notable example of 
use of virtual system can be found in [RS11]. We also remark that Theorem 5 
can be straightforwardly extended to time-dependent switched systems. The 
proof follows exactly the same steps of those used to prove Theorem 5 and 
hence it is omitted here for the sake of brevity. 



Example 1 

As an example illustrating the key features of partial contraction and virtual 
systems, consider a PWSC system of the form 

x = L(t,x)x. (20) 

Notice that such a system may e.g. model a networked control system or a 
network of biochemical reactions. 

We assume that the system is not contracting. That is, the Jacobian matrix 

dL Tl s 
—x + L(t,x), 

does not have any uniformly negative matrix measure. We also assume that 
there exists a uniformly negative matrix measure for L(t,x), i.e. 

3/i : fj, (L(t, x)) < — c, c>0 a.e. 

Clearly, in this case, system (20) is not contracting nevertheless Theorem 5 
can be used to show that, at steady state, all trajectories of (20) converge 
towards a unique solution. In particular, consider the system 

V = v(y,x,t) = L(t,x)y, 

where x, the state variable of the original system, is seen as an external input. 
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It is straightforward to check that 

v(x, x, t) = L(t, x)x = x, 

and hence it is a virtual system in the sense of Theorem 5. Moreover, the 
Jacobian matrix of the virtual system is simply 

dL{t,x)y 
J{t,x) = 7j- = L(t,x). 

Since we assumed that there exists a uniformly negative matrix measure for 
L(t,x), the virtual system is contracting for any x. Therefore, all of its solu- 
tions will converge to a unique trajectory, say y*, such that: 

y* = L(t,x)y* 

Since the solutions of the real system are also particular solutions ofthe virtual 
system, it follows that 

\x(t)-y*(t)\->0 a.e. 

That is, all solutions of the real system will also converge towards y* and, 
hence, towards each other. 



6 Applications 

The extension of contraction and partial contraction to Caratheodory systems 
is a flexible tool that can be used for analysis and design as illustrated by 
means of two representative examples described in this section. 

6. 1 Stability of Piecewise Linear Systems 

Using the concept of contraction for PWS systems, it is possible to easily prove 
the following result to assess the stability of piecewise linear systems of the 
form 

x = A(t,a)x, (21) 

where and a(t) : [0, +oo) — >■ E is the switching signal with S being a 

finite index set. Several stability results for this class of systems are available in 
the literature (see [LA09] for an extensive survey). A classical approach is that 
of finding conditions on the (finite) set of matrices A(t, a) guaranteeing the 
existence of some common quadratic Lyapunov function (CQLF). In [LA09] 
(Theorem 8, p. 311), it is proven that the origin is a globally asymptotically 
stable equilibrium of (21) if and only if there exist a full column rank matrix 
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L E ^rnxn w jth m > n and a family of matrices A4 E jgmxm guc ] 1 
yUoo (Aj) < for all i = a E S. 

Here we show that a related stability condition can be immediately obtained 
by applying contraction theory. Indeed, we can prove the following result. 

Corollary 6 Given a piecewise linear system of the form (21), if the matrices 
A(t, a) are bounded and measurable for any a and there exist some matrix 
measure such that 

n(A(t,a)) < -c, c>0 VtER + , Vcr E S, (22) 

t/ien, all solutions of (21) converge asymptotically towards the origin indepen- 
dently from the switching sequence. 

Proof. Under the hypotheses, system (21) satisfies Theorem 5 and therefore 
is contracting with all of its trajectories converging towards each other. Since, 
x(t) = is also a trajectory of (21), the proof immediately follows. 

As an example, take the system: 

x = A(a)x, a €{1,2} (23) 

with 





(-1.0 






( 


-3.0 


1.0 \ 


A(l) = 






A{2) = 










K 0.8 


-3.0 y 




\ 


2.0 


-l.sj 



Note that using the matrix measure ji\ induced by the 1-norm, we have 
yUi[v4(l)] < and //i[A(2)] < 0. Hence, it is immediate to prove asymptotic 
convergence of all solutions towards each other and onto the origin using Corol- 
lary 6, for arbitrarily switching signal cr(t). 

We wish to emphasize that the result based on contraction embeds as a special 
case the stability condition reported in [LA09]. Indeed, setting L = I and 
Ai = A(a) in Theorem 8, p. 311 in [LA09] is equivalent to using Corollary 6 
with the matrix measure fi^. Moreover, the proof based on contraction can 
also be extended to nonlinear switched systems. 

6.2 Incremental Stability of PWSC and TSS systems 

Results stated in this paper give us a powerful tool to easily show exponential 
incremental stability for switched Caratheodory systems. As an example, we 
consider a switched version of the biological system used in [RdBSlOb] to 
illustrate the application of contraction to biochemical networks. Specifically, 
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following [RdBSlOb], we consider the externally driven transcriptional system 
described by the equations: 



x = u(t) (X tot - x - y) - 8x + fay - fa (e T - y) x, 
y = -fay + fa (e T - y) x. 



(24) 
(25) 



In the above model, x represents the concentration of the transcription factor 
X, y the concentration of the complex protein-promoter Y with the constant 
ex being its total concentration. The external signal u(t) represents the con- 
centration of a transcriptional factor which inactivates transciptor X (and 
consequently X) through a phosphorylation reaction, while the scalars fa and 
fa represent the binding and dissociation rates, respectively. Finally, X to t rep- 
resents the total constant amount of reactant concentrations X tot = X + Y. 

We now assume that the smooth system (24)- (25) is also affected by the 
possible presence of an extra term v(t), of the form: 



From the physical viewpoint, the term v(t) models a degradation (see e.g. 
[SSP06]) on the state variable x which becomes active when the value is above 
a certain threshold, say h. That is, the degradation is activated when the total 
amount of x is greater than h and its effect can be modelled as an additive 
term. 

Using as in [RdBSlOb] the simple change of variables x t = x + y, the resulting 
model can be rewritten as: 



To prove contraction, and hence global incremental stability of this switched 
system, we need to derive the Jacobian which, in this case, is the discontinuous 




x t = u(t) (X tot - x t ) - 8x t + Sy + v(t), 
y = -fay + fa (e T - y) (x t - y) , 



(26) 
(27) 



where the discontinuous term v(t) is given by: 




tfxt~y<h 
/3[x t -y-h] ifx t -y>h 
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function: 



J(x, t) 



J s if x < h 

J s + J ns ifx>h : 



where 



u(t) - 5 



5 



s 



h(e T - y) 



fa + fa(-e T - x t + 2y) 



and 



J, 



ns 







Using /ioo( - ) as a matrix measure, we find that /loo(Js) is negative if the fol- 
lowing inequalities hold: 



for ci, c 2 > 0. 

As shown in [RdBSlOb], the first inequality is always satisfied as the system 
parameters and the periodic input are assumed to be positive. Furthermore, 
taking into account that, for physical reasons, the term e T — y > o, it is easy 
to prove that inequality (29) is also fulfilled. 

We now have to consider the effect of the switching by looking at the mea- 
sure of the matrix J s + J ns . It is immediate to see that, Hoo{J s + Jns) is also 
negative if inequalities (28) and (29) are satisfied. Hence, according to The- 
orem 3, the switched biochemical network under investigation is contracting 
and is therefore incrementally stable. This also implies that, as discussed in 
[RdBSlOb], the transcriptional network continues to be entrainable even when 
an additional discontinuous term is added to the model. This is confirmed by 
the numerical simulation reported in Fig.l. 

Fig.2 shows that, as expected, the switching signals associated to trajectories 
starting from different initial conditions also synchronize asymptotically. 



u(t) - 5+ \S\ < -ci; 



(28) 



fa + k 2 (-e T -x t + 2y) + \k 2 (e T - y)\ < -c 2 ; 



(29) 
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Fig. 1. Convergence of trajectories of the transcriptional module starting from dif- 
ferent initial conditions towards the same unique periodic orbit. Simulation were 
carried out with the following parameter values: k\ = 0.5, &2 = 5, X to t = 1, &r = 1, 
S = 20, P = 1, h = 0.01. The periodic input was set to 1.5 + 2sin(10i). 

6.3 Convergence of networks of time- switching systems 

Contraction analysis can be an invaluable tool to study convergence of net- 
worked systems as proposed in [WS05]. Here, we use the extension of con- 
traction to Caratheodory systems presented in this paper to derive conditions 
guaranteeing convergence of a network of diffusively coupled switched linear 
systems. Specifically, we consider a network of the form: 



where Xi G lR n represents the state vector of node z, Ni denotes the set of 
the neighbors of the z-th node whose cardinality (i.e. the degree of the i- 
th network node) is denoted with <fj. In the above equation, T is a coupling 
matrix, often termed as inner- coupling matrix in the literature. In what follows 
the eigenvalues of the network Laplacian matrix (L) are denoted with A^; A2 
being the smallest non-zero Laplacian eigenvalue (algebraic connectivity). We 
assume A(a(t)) to be bounded and measurable. 

We will now show that, by using contraction, a sufficient condition can be de- 
rived ensuring all the solutions of the network nodes globally exponentially 



A(a(t))xi + r i x j ~ x i\ 



(30) 



jeNi 
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Fig. 2. Synchronization of the switching signals associated to trajectories of the 
modified transcriptional module starting from different initial conditions 

converge, almost everywhere, towards the n- dimensional linear subspace^l 
M. s := {xi = . . . = Xn}. In what follows, we will denote by s(t) the com- 
mon asymptotic behavior of all nodes on M. s . Note that s(t) is obviously a 
solution of each isolated node of (30), i.e. s(t) = A(a(t))s(t). We will also say 
that the network nodes are coordinated (or that the network is coordinated) 
if 

lim \xi(t) - s(t)\ = 0, 



t— >oo 



a.e. 



In the special case where s(t) exhibits an oscillatory behavior, we will say that 
all network nodes are synchronized (or that the network is synchronized). 

Theorem 7 The trajectories of all nodes in the network (30) exponentially 
converge towards each other almost everywhere (i.e., the network is coordinated 
a.e.) if (i) the topology of the network is connected and (ii) there exist some 
matrix measure, \i, such that: 

n(A(a(t))-\ 2 r) <-c, c>0, 

for all a e E and for almost all t. 



1 It is straightforward to check that this subspace is flow invariant for the network 
dynamics. 
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Before presenting the proof of the Theorem, we report here two useful results, 
[Arc]. 

Lemma 8 Let (g) denote the Kronecker product. The following properties hold: 

• (A <g> B) (C <g> D) = (AC) <g> (BD); 

• if A and B are invertible, then (A ® B)~ l = A^ 1 ® 

Lemma 9 For any n x n real symmetric matrix, A, there exist an orthogonal 
n x n matrix, Q, such that 

Q T AQ = U, (31) 
where U is an n x n diagonal matrix. 



Proof (Theorem 7). Define: 



X :-- 



T T 
x 1 , . . . , x N 



S := 1 N (8) s, E := X - S, 



where ljv denotes the iV-dimensional vector consisting of all Is. (Notice that 
such a vector spans M. s .) The network dynamics can then be written as: 

X = (I N ® A(a(t)))X - (L ® T)X, 

so that the error dynamics is described by 

E= (I N ® A{a{t)))E- (L®T)X (32) 

Notice that network coordination is attained if the dynamics of (32) transversal 
to M s is contracting. Furthermore, notice that 

(L ® T)X = (L®T)(E + S) = (L® T)E + (L ® T)S = 

(l ® r)E + (l ® r)(ijv ® s) = (l ® r)F, 

where the last equality follows from Lemma 8 and from the fact that L-l^ = 0, 
since the network is connected by hypothesis. Thus, from (32), we have: 

E= {I N ® A(a(t)))E- (L®T)E. (33) 

Since L is symmetric, by means of Lemma 9 we have that there exist an iV x N 
orthogonal matrix Q (Q T Q = In) such that: 

A = Q T LQ, 

where A is the N x N diagonal matrix having the Laplacian eigenvalues as its 
diagonal elements. 

Now, considering the following coordinate transformation: 

Z=(Q®I n y 1 E, 
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equation (33) can be recast as 

Z ={Q® [(I N ® A(a(t))) - (L <g> r)] (Q ® / n )Z. 



Then, using Lemma 8, we have: 

(Q ® (Jjv <8> (Q ® /„) = 

® 4) (/iv ® ^(*(*))) (Q ® /„) = 
(Q- 1 ®AKt)))(Q®/„) = 

(Jjv ® A(<r(f))) . 

Analogously: 

(Q^Q-^L^T) (Q® /„) = 

(Q- 1 ® /„) (l ® r) (Q ® / n ) = 
(q^l ® r) (Q ® /„) = 

Q- 1 LQ®V = 

a® r. 

That is, the network dynamics can be written as: 

Z=[I N ®A{a{t))-\®T]Z, (34) 

or equivalently: 

Zi = [A(a(t)) - A,r] Zi, i = 1, . . . , N, Zl eR n . 

Now, recall that the eigenvector associated to the smallest eigenvalue of the 
Laplacian matrix, i.e. Ai = 0, is In and spans M. s . Therefore, the dynamics 
along M. s is given by 

il = [A(<7(f))]zi, 

i.e. it is a solution of the uncoupled nodes' dynamics. The dynamics transversal 
to the invariant subspace is given by: 

Zi — [A(a(t)) - XiT] Zi, i = 2,..., N. 

Obviously [A(a(t)) — XiT] is bounded and measurable. Thus, by virtue of 
Corollary 6, all node trajectories globally exponentially converge a.e. towards 
Ai s , if all of the above dynamics are contracting. Now, it is straightforward 
to check that such a condition is fulfilled if 

z 2 = [A(a(t))-X 2 T]z 2 

is contracting. As this is true from the hypotheses, the result is then proved. 
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6.4 A numerical example 



As a representative example, in this Section we use the results presented above 
to synchronize a network of the form (30), where the dynamics of each uncou- 
pled node is given by: 



The matrix T is chosen as: 



Xli 







sin(t)| 


%2i 




-l 









1 




r 


= k 




5 






1 





X j 



(35) 



with k being the coupling gain that will be determined using Theorem 7. The 
network considered here consists of an all to all topology of three nodes. That 
is, 

2 -1 -1 
-1 2 -1 
-1 -1 2 

and A 2 = 3. Thus, from Theorem 7 it follows that the network synchronizes if 
there exist some matrix measure and two positive scalars ci, c 2 such that: 



sin(t) 
-1 

-sin(i) 
-1 





1 


-3k 






1 



3k 



1 
1 



< — Ci, if sin(t) > 0; 
< -c 2 , if sin(t) < 0. 



That is, synchronization is attained if 




< — ci, if sin(i) > 0; 



< -c 2 , if sin(i) < 0. 



Now, using the matrix measure induced by the vector-1 norm (column sums), 
it is straightforward to see that the above conditions are fulfilled if the coupling 
gain is selected as 

*>§■ 
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Fig. 3. State evolution of the network of three switched linear nodes when k = 
(left panel) and k = 0.4 (right panel). 

As shown in Fig. 3, the theoretical predictions are confirmed by the numerical 
simulations. 



7 Conclusions 



An extension of contraction theory was presented to a generic class of switched 
systems: those satisfying conditions for the existence and uniqueness of a 
Caratheodory solution. In particular, it was proven that infinitesimal contrac- 
tion of each mode of a switched system of interest gives a sufficient condition 
for global exponential convergence of trajectories towards each other. This re- 
sult was then used on a set of representative applications. It was shown that, 
by using contraction, it is possible to immediately derive sufficient conditions 
for global stability of switched linear systems. Also, contraction was used to 
obtain sufficient conditions for the convergence of all nodes in a network of 
coupled switched linear systems towards the same synchronous evolution. 

We wish to emphasize that the results presented in this paper can be im- 
mediately applied to generalize to switched Caratheodory systems all of the 
analysis and design results based on contraction analysis available for smooth 
systems in the literature. Examples of applications include nonlinear observer 
design, network protocols design for network coordination, analysis/control of 
asynchronous systems and biochemical systems. 

Finally, note that the results presented in this paper are the first essential 
stage needed to develop a systematic approach to extend contraction analysis 
to generic classes of switched systems. The next step is that of addressing the 
challenging problem of studying convergence in Filippov systems where sliding 
motion is possible. This is currently work in progress and will be presented 
elsewhere. 
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